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the Quark-Meson-Coupling Model 



O 



22 



o 

=5 



X 



J. Rikovska Stone a ' b 

1 Department of Physics, University Oxford, 
Oxford 0X1 3PU, United Kingdom 



O . 

b 'Department of Chemistry and Biochemistry, 
^ ' University of Maryland, College Park, MD 20742, USA 

> ' 

o 



P. A.M. Guichon '* 

c SPhN-DAPNIA, CEA Saclay, F91191 Gif sur Yvette, France 



>\ H.H. Matevosyan d,c 



^Louisiana State University, Department of Physics 
& Astronomy, 202 Nicholson Hall, Tower Dr., LA 70803, USA 



: A. W. Thomas 6 



c Thomas Jefferson National Accelerator Facility, 
12000 Jefferson Ave., Newport News, VA 23606, USA 



Abstract 

A new density dependent effective baryon-baryon interaction has been recently de- 
rived from the quark-meson-coupling (QMC) model, offering impressive results in 
application to finite nuclei and dense baryon matter. This self-consistent, relativis- 
tic quark-level approach is used to construct the Equation of State (EoS) and to 
calculate key properties of high density matter and cold, slowly rotating neutron 
stars. The results include predictions for the maximum mass of neutron star models, 
together with the corresponding radius and central density, as well the properties 
of neutron stars with mass of order 1.4 Mq. The cooling mechanism allowed by 
the QMC EoS is explored and the parameters relevant to slow rotation, namely the 
moment of inertia and the period of rotation investigated. The results of the calcula- 
tion, which are found to be in good agreement with available observational data, are 
compared with the predictions of more traditional EoS, based on the Al8+<5v+UIX* 
and modified Reid soft core potentials, the Skyrme SkM* interaction and a relativis- 
tic mean field (RMF) models for a hybrid stars including quark matter. The QMC 
EoS provides cold neutron star models with maximum mass 1.9-2.1 M©, with cen- 
tral density less than 6 times nuclear saturation density (no = 0.16 fm -3 ) and offers 
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a consistent description of the stellar mass up to this density limit. In contrast with 
other models, QMC predicts no hyperon contribution at densities lower than 3no, 
for matter in /3-equilibrium. At higher densities, and A hyperons are present. 
The absence of lighter S^' hyperons is understood as a consequence of antisym- 
metrisation, together with the implementation of the color hyperfine interaction in 
the response of the quark bag to the nuclear scalar field. 
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1 Introduction 

Nuclear forces play an important role in many stellar environments, acting in 
concert with gravitational forces to form compact objects. For example, the 
observed maximum mass of cold neutron stars cannot be explained without 
considering the pressure, arising from the strong nuclear force, which opposes 
gravitational collapse. It follows that strongly interacting matter is one of the 
fundamental systems that we have to study in order to understand stellar 
phenomena as well as the properties of finite nuclei. 

The properties of stellar matter are strongly density dependent. The theo- 
retical framework for describing the properties of strongly interacting matter 
is the Equation of State (EoS), which relates the pressure and total energy 
density of the system. These quantities are derived from the density depen- 
dence of the energy per particle of the system, calculated using a model for 
the baryon-baryon interaction between particles present in the matter. There 
is a large variety of EoS available in the literature (see e.g. [1,2,3,4]), based 
on very different assumptions concerning the nature of the interactions gov- 
erning the stellar matter. Ways of constraining the available EoS are being 
actively sought (see e.g. [4,5,6,7,8]) but no unambiguous findings have been 
reported as yet. The study of the properties of high density matter and of cold 
neutron stars provides an attractive laboratory for seeking such a constraint. 
However, the fundamental question arises as to the limit of applicability of 
any of the existing EoS at densities as high as 6 4- 10 no, required, for example, 
for calculating the maximum mass of neutron star models. 

At sufficiently high density, nuclear matter will most likely be composed 
not just of nucleons and leptons but also strange baryons and boson con- 
densates [9,10,11,12]. More importantly, it is unlikely that these particles 
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will maintain their identity as confined entities of quarks, but the interac- 
tion between overlapping quark bags will result in partial or full deconfine- 
ment [13,14]. Therefore a realistic high density EoS should contain all this 
physics. Numerical extrapolation of a nucleon-only EoS to high densities, al- 
though it may appear to work (mechanically), is dangerous and may obscure 
some important effects. This is true even at lower densities, in the region 
around twice nuclear saturation density, as some models predict low threshold 
densities for the creation of hyperons in beta equilibrium. The investigation of 
this effect on non-equilibrium matter is expected to be crucial for modelling 
core-collapse supernovae, as the density at bounce after the original collapse 
is around 2 -j- 2.5 no and the presence of hyperons will soften the EoS, thus 
assisting the explosion. 

The relativistic formulation of the QMC model [15,16] offers a unique oppor- 
tunity to self-consistently probe the composition of high density matter. In 
this work, the matter which is investigated is taken to be in /9-equilibrium, 
including the full baryon octet as well as electrons and muons. The density 
and temperature dependent thresholds and production rates of hyperons in 
non-equilibrium stellar matter, relevant for core-collapse supernova physics, 
will be the subject of a forthcoming publication. 

The general QMC model is briefly described in Section 2, followed by a de- 
scription of its application to uniform nuclear matter at zero temperature and 
in /3-equilibrium (Section 3). The basic features of cold, slowly rotating neu- 
tron stars are given in Section 4. Details of QMC neutron star models and 
some comparison with observational data and other model calculations are 
given in Section 5, with a summary of our main conclusions in Section 6. 



2 QMC model 

2. 1 The physics of the model 

The formulation of the QMC model [17] which we use in this work has been 
presented in a previous paper [16]. There it was shown that it is possible to 
build a nuclear Hamiltonian which is consistent with relativity and which can 
be used at high density. The predictions of this model, which represents a 
significant improvement over earlier work [15], because the need to expand 
about < a >= was removed, have been sucessfully confronted with the 
phenomenology of finite nuclei. In this section we rely on the results of Ref. [16] 
and generalize the formalism to nuclear matter containing an arbitrary mixture 
of the octet members, N, A, S, H. 
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We recall that, in the QMC model, the nuclear system is represented as a col- 
lection of quark bags. These bags may contain strange as well as non strange 
quarks, so the framework is well adapted to describing any kind of baryon. 
The crucial hypothesis is that in nuclear matter the baryons retain their in- 
dividuality, or more precisely, they remain the pertinent degrees of freedom 
(or quasi-particles) even in dense nuclear matter. From this point of view the 
usual statement that the bags should not overlap is clearly too restrictive. It 
gives too much weight to the description of the internal structure within the 
bag model, while the latter is in fact used only to infer the interactions be- 
tween the quasi-particles. In this respect, we point out that the bag model is 
an effective realisation of confinement, which must not be taken too literally. 
The QCD lattice simulations of Bissey et al. [18] indicate that the true con- 
finement picture is closer to a Y-shaped color string attached to the quarks. 
Outside this relatively thin string one has the ordinary, non-perturbative vac- 
uum, where the quarks from other hadrons can pass without disturbing the 
structure very much. So, even though the bag model imposes a strict condition, 
which prevents the quarks from travelling through its boundary, this must be 
seen as the average representation of a more complex situation and one should 
not attribute a deep physical meaning to the surface of the cavity nor to its 
size. In particular, estimating the density at which the QMC approximation 
breaks down as the reciprocal of the bag volume is certainly too pessimistic. 
Of course, there is a density above which the QMC model becomes inadequate 
but in this work we assume that this critical density is large enough that we 
can use the model to predict the properties of neutron stars. 

The salient feature of the QMC model is that the interactions are generated by 
the exchange of mesons coupled locally to the quarks. In a literal interpretation 
of the bag model, where only quarks and gluons can live inside the cavity, this 
coupling would be unnatural. Again, we must return to the more realistic 
underlying picture where the quarks are attached to a string but otherwise 
move in the non-perturbative QCD vacuum. There, nothing prevents them 
from feeling the vacuum fluctuations, which we describe by meson fields. As 
before we limit our considerations to the a, u and p mesons. We also estimate, 
in a simplified manner, the effect of single pion exchange. At this point it is 
useful to recall that the a meson used here is not the chiral partner of the pion. 
It is a chiral invariant scalar field, which mocks up the correlated multi-pion 
exchanges in the scalar-isoscalar channel. 



2.2 Effective mass and energy of a bag in the nuclear medium 

In Ref. [19] we used the Born-Oppenheimer (BO) approximation to derive 
the energy of a quark bag coupled to the nuclear mean fields associated with 
the ex, us and p mesons. Given the position and velocity of a bag, we obtained 
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Table 1 

Quantum numbers of the octet members. 

its energy by solving the field equations for the confined quarks. As the BO 
approximation requires that both position and velocity are known, this energy 
is of course classical and canonical quantization is performed after the full 
Hamiltonian has been built. As our previous work was limited to nuclear 
matter and finite nuclei built primarily of nucleons, we now generalize our 
results to include hyperons. Here we limit our considerations to the spin 1/2 
SU(3) octet (N, A, S,S) and therefore a baryon can be specified by \f >= 
\tms >, with t,m the isospin and its projection and s the strangeness - see 
Table 1. Our working hypothesis is that the strange quark is not coupled to the 
meson fields, which can be seen as a consequence of the fact that the mesons 
represent correlated pion exchanges. This is also the natural explanation for 
the small spin-orbit splitting observed in A hypernuclei, given that all the A 
spin is carried by the strange quark. We set the u and d quark masses to zero 
and assume that their coupling to the meson fields respects isospin symmetry. 
Then the energy of a baryon of flavour / with position R and momentum P 
has the form|_U 

E = y JP* + M f (ay + glu + g p b.I i (1) 

where Mf(a) is the effective mass, that is the energy of the corresponding 
quark bag in its rest frame, in the a field, and P is the isospin operator for 
isospin t, defined as 

< tms\ J2 T i\t'm's' >= 5(tt')5(ss')il m) (2) 



with T q the Pauli matrix acting on the u and d quarks. Apart from the p 
coupling, for which we denote the field b a , with a = (—1,0,1) the isospin 
index, the energy is diagonal in flavour space. Note that the flavour dependence 
of the p coupling is entirely contained in /*. In Eq. (1) the meson fields are 
evaluated at r — R, the center of the bag0 and, in comparison with Ref. [16], 
we have dropped the spin-orbit interaction - since in the following we shall 

1 We use the system of units such that h = c = 1. 

2 The variation of the field over the confinement cavity generates the spin-orbit 
interaction, as derived in Ref. [16]. 
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consider only the case of uniform matter. The flavour dependence of the oj 
coupling is fixed by the number of non-strange quarks in the baryon, that is: 



£ = uffQu = (1 + - (3) 
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with g u the u — N coupling constant. As shown in the Appendix, the effective 
mass is well approximated by a quadratic expansion 

M f (a) = M f - w° f g a a + ^w° f (g a a) 2 (4) 



where g a is the free a — N coupling constant. The scalar polarisability, d, de- 
scribes the response of the nucleon to the applied scalar field. It is at the origin 
of the many-body forces in the QMC model. The corresponding scalar polar- 
isability for baryon / is dwZ(wZ)~ 2 . The weights wZ,wZ control the flavour 
dependence and in first approximation wZ ~ wZ ~ l + s/3. The hyperfine 
color interaction breaks this relation and the exact values, which depend on 
the free nucleon radius, R{^ ee , are given in the Appendix. 



2.3 Hamiltonian 



The total energy of the nuclear system is the sum of the baryon energies (1) 
and of the energy stored in the meson fields: 



tot 



E 



mesons 



i=l,A 

df 



Er, 



{Vaf + m 2 a a 2 



-- I df 



-Idf 



2i72 

m p b 



(5) 



with m a , m^, m p the meson masses. As usual, we consider only the time compo- 
nent of the vector fields. By hypothesis, the meson fields are time independent. 
Therefore the classical Hamiltonian of the nuclear system is simply 



H(Ri, Pi) = E tot (Ri, Pi, a -> a soU uj — >■ u so i, b —> b 



sol J 



(6) 



where a so i, u so i, b so i are the the solutions of the meson equations of motion: 

5E _5E _ 5E 
5a 5uj 5b n 
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The equations for the u, b fields, which are linear, present no difficulty. By 
contrast, the a field equation is highly non-linear, because of the scalar po- 
larisability term in the effective mass. In Ref. [16], we proposed to solve it 
approximately by writing: 

a = a + 5a, a =< a > (8) 



where, < • ■ • >, denotes the nuclear ground state expectation value, and we 
considered the deviation, 5a, as a small quantity. We refer to Ref. [16] for the 
details of the derivation and here we simply recall the results. The piece of 
the Hamiltonian which depends on a may be written: 



df 



K(a) 




(9) 



where normal ordering is implicit. The one body kinetic energy operator, K(a), 
is defined as (V is the normalization volume): 

k,k'j 



with agy the destruction operator for a baryon of momentum k and flavour /, 
the spin label being omitted for simplicity. Note that the mean field approxi- 
mation amounts to setting 5a = in the Hamiltonian (9). 

The meson field solution in the general case is of little interest to us, so we 
directly give the solution for a uniform system in which case a is a constant 
determined by the self consistent equation: 

1 dK 

m z a oa 



(11) 



which must be solved numerically and the fluctuation 5a is given by: 



where the effective a mass is: 

ml = ml+< — >. (13) 
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2.4 Energy density in the Hartree Fock approximation 



The last step in our formal development is to use Eqs. (9-12) to compute 
the energy density of nuclear matter in the Hartree-Fock approximation. The 
ground state of the system is specified by a set of Fermi levels, {kp(f),f = 
p,n...}, and their corresponding baryonic densities {%} with rig = J2f n f- 
One has: 

M/) 

<*(*)>= (2^)3 E / My/k* + M}{a) (14) 



and the energy density takes the form 




(2tt^ 



M/) 

J dkidk 2 



f o 



kf + M^a) — Jk$ + M*(a). (15) 



(h - k 2 ) 2 +rhldo-v da 



Finally we must add the contributions < V w > and < V p >, corresponding to 
uj and p exchange. The derivation follows the same line as for the a exchange 
but is much simpler and gives the exact solution because these interactions 
are purely 2-body . One gets: 



k F (f) 2 

■GTklV -tt; I dhdk 2 ^ ^ (16) 



<V P > 



V 



^ \tms I 



kp(tms) kp(tm's) 

r -> r m z r 



tram' 



-Gp E ^-iLnJ^E J ^ J ^^-k^ + ml 



As usual we define: 



The a, lo and p mesons represent the multi-pion exchanges that we believe to 
be the most relevant ones. However, this cannot take into account the long 
ranged single pion exchange. In the Hartree Fock approximation, the latter 
contributes only through the exchange (Fock) term, so its contribution to 
the energy is modest. Nevertheless, we think it is worthwhile to study its 
impact on our results because, as a consequence of its long range, the pion 
contribution to the energy has, as compared to the heavy mesons, a specific 
density dependence. With this in mind, we generalize the expression of Ref. [20] 
so as to include the contribution of the full octet. We get: 



<V n > 



9a 



V n B \2f^ / 
16 



24 

J pp + 4 J pn + J nn — — ( Ja,£- + <^A,E° + ^A,£M 



+ 2^ (^s-s- + 2J E - E o + 2J E + s o + Je+^ 
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4J= 



./-()-!) 



(19) 



with : 



M/)M/') 

In Eq. (19), qa = 1-26 is the axial coupling constant of the nucleon, f n = 
93 MeV is the pion decay constant and m n its mass. In the loop integral 
Jff, the first term is the so-called contact term. In practice it can hardly be 
distinguished from the short-ranged contributions of the heavy mesons and 
we omit it. We have checked that it would only induce a small readjustment 
of G a , Gu, G p . In the same spirit we have omitted any 7r-baryon form factor in 
Jff, because the effect of a form factor is to cut off the large momenta of the 
loop and it therefore amounts to a short ranged contribution. In summary, we 
keep only the long ranged Yukawa piece of the pion exchange. 



(p — p 



mi 



(20) 



2.5 Fixing the parameters 



The parameters of the model are the couplings G a , G u , G p , the meson masses 
and the free nucleon radius. We set the n, to, p masses to their physical values. 
We have tried several values of the free nucleon radius, R^ ee and found that 
this has very little influence on our results. We therefore take R^ ee = 0.8 fm, 
which is a realistic value [21]. The a meson mass is not well known because 
the width of the ttti resonance is very large in the physical region. In our study 
of finite nuclei [16], we found that m a ~ 700 MeV produced the best results, 
so this will be our prefered value. However, the sensitivity of the calculations 
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of finite nuclei to the a meson mass comes from the fact that it controls the 
shape of the nuclear surface, a factor which is irrelevant in the case of neutron 
stars. Moreover, since the calculations for finite nuclei are based on a non- 
relativistic, low density approximat iorf 3 ! to the Hamiltonian (9), we must be 
conservative about the a meson mass when we use it in the high density region. 
Therefore, to quantify the influence of this mass, we shall show some results 
for m a = 600 MeV. 

The remaining free parameters, G a , G w and G p are adjusted to reproduce 
the binding energy and the asymmetry energy of ordinary nuclear matter at 
the saturation point. Since we want to sample the stiffness of our equation 
of state, we allow some freedom in the way the couplings are fixed. This is 
summarized in Table 2. First we set the pion contribution to zero and adjust 
the couplings to reproduce the experimental values of the binding energy, (8 = 
— 15.865 MeV), and asymmetry energy, (a s = 30 MeV), of nuclear matter 
at the saturation point (no = 0.16 fm -3 ). This defines the models QMC600 
and QMC700, according to the value of the a meson mass. For all the other 
models (QMCVz, i — 1,4) we set m a = 700 MeV and the effect of the pion is 
included. As can be seen in the last column of Table 2, the incompressibility 
modulus of both QMC600 and QMC700 is above the accepted range, 200 
300 MeV. The most efficient way to decrease is to add an attractive 
interaction which has a weak density dependence at the saturation point. 
The pion Fock term, which is attractive and behaves roughly as p 1//6 , is thus 
a natural candidate. This is confirmed by the model QMC7rl, where the pion 
contribution is included and where drops from 340 MeV to 322 MeV. A 
trivial way to mock up the effect of an almost constant attractive interaction 
is simply to decrease the absolute value of the experimental binding energy. 
By setting £ = —14.5 MeV, a rather moderate distortion of the experimental 
number, we already find = 301 MeV. This defines QMC7r2. In order to 
sample a wider range of values of the stiffness, we allow ourself to artificially 
vary both the pion contribution and the binding energy. This leads to model 
QMC7r3 (resp. QMC7r4), where the pion contribution is multiplied by 1.5 
(resp. 2) and the binding energy is set to —14 MeV (resp —13 MeV). This 
gives Koo = 283 MeV and = 256 MeV, respectively. Note that increasing 
the pion contribution by a factor of 2 is not so arbitrary, because we know 
that it is amplified by RPA correlations [22]. On the other hand, changing the 
experimental binding energy is more difficult to justify and we take it only as 
a crude way to sample the stiffness of the model, given its eventual application 
in the high density region, far away from the point where S is measured. 



3 By this we mean a density about or smaller than the density of ordinary nuclear 
matter. 
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Table 2 

The couplings for different versions of the model. The column ir is the number 
by which the pion contribution has been multiplied. £ is the binding energy of 
symmetric nuclear matter and its incompressibility modulus. 

3 Cold Uniform Matter 



In this section we determine the equation of state of uniform matter, which 
is directly relevant for studies of the interior of neutron stars. The structure 
of the neutron stars is discussed in the next section. Uniform matter in cold 
neutron stars is in a generalized beta-equilibrium, achieved over an extended 
period of time. All baryons of the octet can be populated by succcessive weak 
interactions, regardless of their strangeness [23] . We consider matter formed by 
baryons of the octet, electrons and negative muons, with respective densities 
{n f ,f =p,n, ...} and {n e ,n^}. 

The equilibrium state minimises the total energy density, e, under the con- 
straint of baryon number conservation and electric neutrality. We write: 

e = e B (n p , ...) + e e (n e ) + e^(n^) (21) 



where the baryonic contribution 

e B {n p , ...) = (22) 



is calculated from Eqs. (15-19). It is related to the binding energy per baryon, 
£, by: 

e B (n p ,...)=Y.( £ + M f) n f- ( 23 ) 
/ 



The baryonic pressure, (Pb), the incompressibility modulus, (i^oo), and the 
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sound velocity, (y s ), of baryonic matter have the following expressions: 



2 9 e B ^ n dP B j n B K 



Pb = n B - ± = 9-^, v s = J - (24) 

dn B n B dn B V 9(e B + P B ) 



where n B = J2f n f an d the derivative, with respect to n B , is taken at constant 
fractions (rif/n B , f = 1, 8). Similar expressions hold for each lepton. In proton- 
neutron matter another important variable is the symmetry energy, S, defined 
as the difference between pure neutron and symmetric matter : 

S(n B ) = £{{n p = 0, n n = n B ) -S(n p = n B /2, n n = n B /2) . (25) 



For the energy density of the lepton /, of mass mi and density n;, we use the 
free Fermi gas expression: 

- M,) 4(0 



ej(ni) = -j—-^ J dkyjk 2 + mf, withn* = (26) 



The equilibrium condition for a neutral system with baryon density n B is 



5[e B {n p , ...) + e e {n e ) + e M (n M ) 

71 f ~ u b) + "(£ n f q f - (n e + n M )] = (27) 
/ / 

where qf is the charge of the flavor /. The constraints are implemented through 
the Lagrange multipliers (A, v) so the variation in Eq. (27) amounts to inde- 
pendent variations of the densities, together with variations of A and v. If one 
defines the chemical potentials as 

=w,< =k - (28) 



the equilibrium equations become: 



Hi + A + vq { : 
[i e -v- 



I 



f 



n. 



0. 
0. 
0. 

0. 
0. 



f = p, 



(29) 
(30) 
(31) 
(32) 

(33) 
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This is a system of non-linear equations for {n p , ...,n e , n M , A, u}. It is usual to 
eliminate the Lagrange multipliers, (A, u), using Eqs. (30) and (29), with / = 
neutron. However, for a given value of rig, the equilibrium equation, Eq. (27), 
generally implies that some of the densities vanish and therefore that the 
equations generated by their variation drop from the system (29-31) - because 
there is nothing to vary! In particular, substituting v = fi e in Eqs. (29) is not 
valid when the electron disappears from the system. The equations obtained 
by this substitution may have no solution in the deleptonized region, since one 
has forced v = 0. To correct for this, one must restore v as an independent 
variable when one reaches this region. This is technically inconvenient and we 
found it is much simpler to blindly solve the full system of equations, (29-33), 
for the set {n p , n e , n^, A, v}. The only simplification which is not dangerous 
is the elimination of the muon density in favor of the electron density by 
combining Eqs. (30,31) to write /x^ = /x e , which is solved by 

k F (/i) = &yfk F (e) 2 + ml-ml (34) 



where 3ft denotes the real part. This is always correct because if the electron 
density vanishes then so does the muon density and the relation (34) reduces 
to = 0. 

To solve the system (29-33), let us define the set of relative concentrations 
(note that the lepton concentrations are also defined with respect to 

Vn B n B n B n B n B J 



Each member, y$, of the set Y is associated with an equation Ei , the one among 
Eqs. (29,31) which is obtained by taking the variation of e with respect to pi = 
Uin B . Let us assume that a solution, Y , has been found at some baryon density, 
n B = no. One first tests whether the threshold for some species, i, is crossed 
when uq is incremented by 8n. Since yi — below and across the theshold, the 
condition for the appearance of the species is that Ei(n B , Y ) changes its sign 
between n and n +5n. If this happens, the equation Ei is added to the system. 
If the threshold condition is met for several species, all the corresponding 
equations are added. The system is then solved numerically at n B = no + Sno, 
using Yq as a first approximation. When a concentration drops below a given 
small number, n, it is set to zero and the corresponding equation is removed 
from the system. The value of n depends on the accuracy of the solution and 
in our calculation it has been set to 10 -4 . We have checked that rj = 10~ 3 
gives the same result. We solve the system with the initial condition that 
at n B = the matter contains only neutrons. Once the equilibrium solution 
Y{n B ) has been found for the desired range of baryon density, typically n B = 
0-rl. 2fm~ 3 , it is used to compute the equilibrium total energy density, e(n B ) = 
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^{VpfiBi UuMb)- The corresponding total pressure, P(n B ), is computed as 
the sum of the baryon pressure, (24), and of the lepton pressures. Using the 
equilibrium equations, Eqs. (29-33), it is straightforward to show that the total 
pressure can also be evaluated as (note the total derivative): 

P[n B ) = n B - . 36 

dn B n B 



As a numerical check, we used both expressions, which did indeed yield the 
same result within the numerical accuracy. 



4 Neutron Stars 



A neutron star is an object composed of matter at densities ranging from that 
of terrestrial iron up to several times that of nuclear matter. In cold neutron 
st cirs, £is the baryon density tib increases from about 0.75 no up to 2 -r- 3no, 
stellar matter becomes a homogeneous system of unbound neutrons, protons, 
electrons and muons and, if enough time is allowed, will develop /3-equilibrium 
with respect to the weak interactions. All components that are present on a 
timescale longer than the life-time of the system take part in equilibrium. For 
example, neutrinos created in weak processes in a cold neutron star do not 
contribute to the equilibrium conditions as they escape rapidly. At even higher 
densities, heavier mesons and strange baryons play a role (see e.g. Ref. [24] 
and references therein, [25,26,27,28]). Ultimately, at the center of the star, 
a quark matter phase may appear, either alone or coexisting with hadronic 
matter [23,29,30,31]. Below ub ~ 0.75 no, the matter forms the inner and 
outer crust of the star, with inhomogeneities consisting of nucleons arranged 
on a lattice, as well as neutron and electron gases. In this density region the 
QMC Equation of State (EoS) needs to be matched onto other equations 
of state, reflecting the composition of matter at those densities. The Baym- 
Bethe-Pethick (BBP) [32] and Baym-Pethick-Sutherland (BPS) EoS [33] are 
used in this work. 

Setting up an EoS over the full density range, up to at least 6 n , allows 
calculation of one of the most important observables of neutron stars, the 
maximum gravitational mass, M™ ax , and the corresponding radius, i? m ax- The 
most accurately measured masses of neutron stars were, until very recently, 
consistent with the range 1.26 to 1.45 M Q [5]. However, Nice et al. [34] recently 
provided a dramatic result for the gravitational mass of the PSR J0751+1807 
millisecond pulsar (in a binary system with a helium white dwarf), M g = 2.1± 
0.2 M (with 68% confidence) or 2.li^5 with 95% confidence which makes it 
the most massive pulsar measured. This observation potentially offers one of 
the most stringent tests for the EoS used in the calculation of cold neutron 
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stars. It also sets an upper limit to the mass density, or equivalently, the energy 
density, inside the star [5]. A lower limit to the mass density can be derived 
using the latest data on the largest observed redshift from a neutron star, 
combined with its observational gravitational mass. 

The calculated maximum mass is determined to large extent by the high den- 
sity EoS. It has been argued that extrapolation of the nucleon-based EOS, 
built using an effective force (e.g. non-relativistic Skyrme [4], relativistic mean 
field [23]) or phenomenological interactions (such as Al8+5v+ UIX* (APR) [35]) 
to densities corresponding to M™ ax , may not be unreasonable [35,36] and that 
the error made by not including the heavy baryons and possible quarks in the 
calculation may not be significant. On the other hand, we note that in the 
hybrid calculations of Lawley et al. [13], the maximum mass of a neutron star 
was significantly reduced. In any case, we will explore the consequences of this 
extrapolation by using the QMC EoS and taking into account the different 
composition of the homogeneous matter, starting from a nucleon-only model 
and then including the full baryon octet. We shall also comment shortly on the 
presence of exotics, such as penta-quarks and 6-quarks bags, in high density 
matter. Furthermore, neutron star models at around the 'canonical' mass of 
1.4M , with central densities of the order 2-=-2.5no and lower, are also studied 
as representative cases for nucleon-only based EoS. 



4-1 Cold non-rotational neutron star models 



The gravitational mass and the radius are calculated using a tabulated form of 
the composite EoS with a chosen QMC interaction. The Tolman-Oppenheimer- 
Volkov equation [37,38] 

dP _ Gmp (1 + P/p) (1 + 4nr 3 P/m) 
dr r 2 1 — IGmjr 

is integrated with 

m{r) = J Airr' 2 p(r') dr', (38) 
o 



where p is the total mass density and G is the gravitational constant, in or- 
der to obtain sequences of neutron-star models corresponding to any specified 
central mass density. The solution gives directly the radius, R, of the star (the 
surface being at the location where the pressure vanishes) and the correspond- 
ing value for the total gravitational mass M g = m(R). 
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It is also important to calculate some other properties of these neutron star 
models. The total baryon number A is given by 

R 

4:nr 2 riB{r)dr 



A = / avJ 1/2 . (39) 

1 - 2Gm{r)/r) 1/2 



The total baryon number, A, multiplied by the atomic mass unit, 931.50 MeV, 
defines the baryonic mass M . The binding energy released in a supernova 
core-collapse, forming eventually the neutron star, is approximately 

E bmd = (Am - M g ) , (40) 



where mo is defined as the mass per baryon of 56 Fe. Analysis of data from 
supernova 1987A leads to an estimate of Eund = 3.8 ± 1.2 x 10 53 erg [39]. 

The gravitational red shift of the photons emitted radially outwards from a 
neutron star surface is given by 

Zsurf = ( 1 - — " 1 (41) 



It is an important quantity needed, for example, for obtaining the mass and 
radius of a neutron star separately [5]. Other quantities of interest for pos- 
sible comparison with observational data are the minimum rotation period, 
Tmin [40], and the moment of inertia X (see Ref. [6] and references therein). The 
minimum period is given by the centrifugal balance condition for an equatorial 
fluid element (i.e., the condition for it to be moving on a circular geodesic). 
While determining this accurately requires using a numerical code for con- 
structing general-relativistic models of rapidly rotating stars, quite good val- 
ues can be obtained from results for non-rotating models using the empirical 
formula [41,42] 

/ /If max \ -1/2 . td x 3/2 

r mm = 0.820 ( ) ms . (42) 

V M Q J VIOkm/ v ; 



The shortest period so far observed is 0.877 ms [40] but it is possible that 
this limit may be connected with the techniques used for measuring pulsar 
periods, rather than being a genuine physical limit. 

The moment of inertia X is calculated in this work by assuming that the 
rotation is sufficiently slow that the general relativistic slow rotation approxi- 
mation [43] gives a rather accurate description of the star and the surrounding 
space-time and using the expressions given in Ref. [44]. 
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5 Results and discussion 



5.1 Uniform matter 

The properties of uniform matter and neutrons stars were calculated for 12 
EoS designed to systematically investigate the dependence on the basic QMC 
model parameters, summarised in Table 2. There are three types of EoS used 
in this work. They have the same QMC parameters but differ by the limitation 
that we impose on the matter composition. The sets F-QMCx, where x refers 
to the parametrisations of Table 2, describe matter consisting of the full baryon 
octet plus electrons and muons. The sets N-QMCx are used for n + p + e + \x 
matter and PNM-QMCx for pure neutron matter. For simplicity, we only 
show results for QMC700 and QMC7r4 in this section. These sets have the 
most realistic values of the QMC parameters and correspond to two rather 
different values of 340 and 256 MeV. In the next section all the QMC 
parametrisations of Table 2 are used for the calculation of neutron star models 
and their properties are discussed in more detail. 

In addition to the QMC EoS, we use four other models for a comparison 
based on very different physics. The older Bethe- Johnson (BJ) EoS [45] is 
based on the modified Reid potential and the medium effects are treated us- 
ing constrained variational principle (model C in [24]). The BJ EoS describes 
uniform matter including n,p, A, E ±,C) , A ±,Q and A ++ . The more recent EoS, 
calculated with the Al8+5w+UIX* potential (APR) [35], is based on the Ar- 
gonne A18 two-body interaction and includes three-body effects through the 
Urbana UIX* potential. It also includes boost corrections to the two-nucleon 
interaction which corresponds to the leading relativistic effects at order v 2 . 
This interaction is considered one of the most modern phenomenological po- 
tentials used for the description of nuclear matter containing nucleons only. 
As an example of a phenomenological, nucleon-only EoS the Skyrme SkM* 
interaction is used, in particular because of the similarity of its results for fi- 
nite nuclei to those of the effective force derived from the QMC model [15,16]. 
Finally, we also compare with the EoS of a hybrid (neutron-quark) star matter 
of Glendenning (Table 9.1 of [23]) (Hybrid) calculated in RMF approximation. 
The binding energy per particle in pure neutron matter (PNM) and symmet- 
ric nuclear matter (SNM), for QMC700 and QMC7r4, are illustrated in Fig. 1 
in comparison with the APR and SkM* predictions. At lower densities the 
calculated values are very similar. However, with increasing density the differ- 
ences become more dramatic. In particular, while the QMC and APR models 
predict a steady increase of the energy per particle in both pure neutron and 
symmetric matter at about the same rate, the SkM* model yields a more rapid 
growth of the energy per particle in SNM than in PNM, leading to crossing 
at about 0.9 fm~ 3 . As shown in Fig. 2 and discussed in detail in [4], this leads 
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Fig. 1. The energy per particle for pure neutron and symmetric nuclear matter is 
plotted as a function of baryon number density hb for the QMC700 and QMC7r4 
interactions (left panel). For comparison, also shown are equivalent curves calculated 
using the APR potential [35] and the Skyrme SkM* parametrisations (right panel) 



to negative symmetry energy and a collapse of nuclear matter. Although the 
QMC and Skyrme energy functionals have similar structure and yield similar 
results for finite nuclei at normal nuclear density, QMC has a different density 
and isospin dependence [16] that improves the EoS at high densities in com- 
parison with SkM*, while the symmetry energy shows a slower increase with 
density than the APR model. There is no experimental evidence concerning 
the density dependence of the symmetry energy, as it is not possible to study 
the properties of uniform matter at high densities in terrestrial conditions. 
Even in heavy-ion collisions at high energies the maximum baryon density 
reached is unlikely to exceed 2 3no in the phase of the reaction closest to the 
conditions in cold neutron stars. Nevertheless, there is some indirect evidence 
which seems to support models which predict that the symmetry energy grows 
with density [4] . Such models yield sensible results for properties of neutron 
stars that can be compared with observations (see Section 4). 

The results of our calculation of the equilibrium composition of uniform matter 
for QMC700 and QMC7r4 are shown in Fig. 3. In a striking difference from 
all other models which include the full baryon octet, QMC does not predict 
the production of hyperons at densities less than 3 times nuclear saturation 
density. Furthermore, the first hyperons to appear are the cascades, S, together 
with the A hyperon. The S hyperons are not produced at densities below 1.2 
fm -3 . This scenario is a direct consequence of three factors which are present 
in the QMC model and absent in the others. First, we use the Hartree-Fock 
approximation, while the earlier calculations worked only at the Hartree level. 
Actually, we see no reason to ignore antisymmetrisation, which is the key to 
our successful interpretation ofthe spectra of finite nuclei [16]. Second, the 
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Fig. 2. Symmetry energy as a function of baryon number density ns for QMC700 
and QMC7r4. Also shown are the corresponding results for the APR potential and 
the Skyrme SkM* parametrisations. 
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Fig. 3. Relative population of baryons and leptons in /^-equilibrium matter for 
QMC700 and QMC-7r4. For more explanation see text. 

color hyperfine interaction, which is necessary to explain the mass splitting 
within the octet, as well as the NA splitting, induces a flavor dependent 
effective mass which strongly disfavors the E. Third, the scalar polarisability, 
d, induces many body effects which tend to amplify the previous ones. 

To make these considerations quantitative, we have computed the single parti- 
cle threshold densities on top of neutron matter, using the QMC700 parametriza- 
tion and the variations on it defined earlier. The results are shown in Table 3, 
where the first column corresponds to the full QMC700 model as shown in 
Fig. 3. In the second column we show the same calculation but with the hy- 
perfine color interaction and the scalar polarisability set to zero for the strange 
particles. Clearly this produces a very spectacular effect on the £ thresholds, 
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Table 3 

Single particle threshold densities (in fm -3 ) on top of neutron matter for QMC700, 
either in Hartree or Hartree-Fock approximation and with the hyperfine color in- 
teraction (hfc) and scalar polar isability, d, set to zero in the strange sector. 

which become comparable to the S thresholds. In the third column we show 
the same calculation as in the first but now in the Hartree approximation, that 
is dropping the exchange terms in Eqs. (15-17) and readjusting the coupling to 
reproduce the saturation point of nuclear matter. One observes again a strong 
rearrangement of the threshholds in favor of the S. Finally, the last column 
shows the combination of both effects. With these three alterations, the QMC 
model obviously becomes similar to the SU(3) relativistic mean field models 
- with respect to strangeness production. In particular, the first threshold is 
now for the S~, almost identical to that for the A and it occurs at less than 
twice the nuclear density. 



As expected, the presence of hyperons does soften the EoS, as demonstrated 
in Fig. 4. We observe that QMC700 and QMCV4 both show a very similar 
EoS to that for the BJ model, which also contains hyperons. However, it is 
clearly demonstrated that the composition of the stellar matter is not the only 
way to soften the EoS. The pressure also depends critically on the potentials 
acting between the particles present and, in turn, on the density dependence 
of the symmetry energy. Thus, for example, the SkM* Skyrme model pro- 
duces a very soft EoS, in contrast to APR, and almost the same EoS as the 
neutron+quark RMF hybrid model. Clearly more constraints are needed to 
distinguish amongst these possibilities. Further illustrations of the variety of 
possibilities for the dependence of the pressure on the total energy density can 
be found in Ref. [46]. 

To conclude this Section, let us mention that we have also studied the possi- 
bility that matter at high density may contain exotic particles, such as pen- 
taquarks or six quarks bags. The calculation goes along the same lines as for 
the ordinary baryonic matter and poses no conceptual difficulties because, in 
beta equilibrium, one does not worry about the exact way a particular compo- 



20 





1 * 1 




No 


F2, 


r 


(Z 




n 

u 




o.uzo iu 








N-QMCvr4 







4.634 10~ 3 


1.623 


346.7 


192.1 


N-QMC700 


2.62 10 - 7 


3.197 


0.0251 


1.286 


522.1 
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Table 4 

The parameters of the fit defined by eq.43 

nent has been produced. The practical problem is that these would-be exotics 
have not been identified as free particles and therefore their free masses (if 
they exist) are essentially unknown. To get an estimate of these masses, which 
are of course the critical parameters, we must rely on the bag model. It turns 
out that, even if we allow the exotics to be lighter than the bag model predic- 
tion by as much as 400 MeV, they do not appear at any significant level when 
the conditions for beta equilibrium are imposed on matter where the full octet 
is allowed. 



For the readers convenience we have preformed an analytical fit of our EOS 
for the cases shown on Figure 4. It has the form: 



P 



N^ 1 



N 2 e 



1 + exp[(e — r)/a] 1 + exp[— (e — r)/a] 



(43) 



The parameters (Ni,pi, N 2 ,p 2 , r, a) are given in Table 4 with (P, e) expressed 
in MeV I fm 3 . The fit is accurate in the range ps = 4- 1.2/m -3 . 



5.2 Neutron Stars 



The calculated properties of neutron stars with the maximum mass for all 12 
model EoS are summarised in Table 5. For comparison, the results for the 
BJ, APR, SkM* and Hybrid EoS are also included. The same information is 
given for these models for stars with 'canonical' mass 1.4 M in Table 6. It 
should be noted that the entries in Table 6 do not differ for the full octet 
and nucleon-only versions of the model, because the central density of 1.4 M Q 
stars is below the threshold for the creation of hyperons. 

The relation between the gravitational mass and radius for various neutron 
star models is shown in Fig. 5, for a representative set of EoS, chosen to be 
F-QMC700, F-QMCvr4, N-QMC700 and N-QMCvr4. 

Apart from the basic mass-radius relation, there are other features of neutron 
stars that provide the possibility to constrain the EoS and thus the nucleon- 
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Fig. 4. Left panel: Pressure as a function of total energy density for sample EoS for 
/^-equilibrium matter containing nucleons p+n+e+/x (N-QMC700, N-QMC7r4) and 
the full octet (F-QMC700, F-QMCvr4). Right panel: Results of BJ, APR, SkM* and 
hybrid EoS. The curve P = e shows the stiffest possible EoS. 
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Fig. 5. The gravitational masses of non-rotating neutron-star models (measured in 
solar masses) plotted against radius (in kilometers), calculated for selected QMC 
EoS - see the text for more explanation. 

nucleon interaction in stellar matter. The volume and quality of the obser- 
vational data has increased in recent years, offering more options for testing 
the physical basis of the EoS. An important constraint on the EoS at rel- 
atively low densities has been recently identified [7] in connection with the 
very precise determination of the gravitational mass of Pulsar B in the system 
J0737-3039, M 6 = 1.249±0.001 M . If the progenitor star was a white dwarf 
with O-Ne-M core and Pulsar B was formed in an electron-capture supernova, 
a possibility supported in part by new observations [48] , a rather narrow range 
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Table 5 

Parameters of maximum mass, cold neutron stars for the full set of QMC EoS. n c 
and p c are the central baryon number density [ fm~ 3 ] and total mass density [10 15 g 
cm -3 ] respectively, R [km] is the radius, the maximum gravitational mass is M©, A 
[10 57 ] is the total baryon number, E^ind [10 53 erg] is the binding energy and r [msec] 
is the estimate of the rotation period calculated from Eq. (42), v s is the speed of 
sound and z sur f is the radiational red-shift. The index 'max' has been dropped for 
all observables in the header of this table. 
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for the baryonic mass, Mo, between 1.366-1.375 M©, can be determined for 
the pre-collapse core. Assuming that there is no (or negligible) baryon loss 
during the collapse, the newly born neutron star will have the same baryonic 
mass as the progenitor star. For any given EoS for cold neutron star matter 
the relation between the gravitational and baryonic mass can be calculated 
and tested against a very narrow window (full black line box in Fig. 6)de- 
fined by the data for Pulsar B. The width of the window reflects uncertainty 
in modelling the composition of the progenitor core, namely electron-capture 
rates, nuclear network calculation, Coulomb and general relativity corections. 
Very recently, another simulation of the explosion of O-Ne-Mg cores [47] has 
been performed which predicts the baryonic mass of Pulsar B 1.36±0.002 M© 
(full blue line box in Fig. 6). The error includes the uncertainty in the EoS 
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Table 6 

The same as Table 5 but for 1.4 M star models. Data for N-QMCx EoS are not 
shown as they are identical to results of F-QMCx EoS because the central density 
of all 1.4 Mq stars is predicted to be below the threshold for the appearance of 
hyperons. 
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Fig. 6. Relation between the gravitational mass, M g , for various neutron star models 
(including the QMC EoS) and the corresponding baryonic mass, Mo- The boxes rep- 
resent constraints derived by Podsiadlowski et al.[7] (black box) and more recently 
by Kitaura et al.[47] (blue box), based on the properties os system J0737-3039, as 
discussed in the text. 

used and the wind ablation after the simulation was terminated. This model 
(see Ref. [47]) includes a small mass loss between 0.014 to 0.017 M & in the 
process. As seen in Fig. 6, F-QMC7r4 (the other versions give essentially the 
same curves) agrees with the latter constraint and predicts mass loss in the 
region of 0.008 - 0.017 M Q in the model of Podsiadlowski et al [7], which is 
rather similar to the result of Kitaura et al [47]. The other sample EoSs either 
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Fig. 7. Moment of inertia, scaled by M 3 / 2 , as a function of the stellar mass, M, 
for sample QMC EoS. The BJ, APR, SkM* and Hybrid EoS results are added for 
comparison. The dash-dot lines indicate the predicted experimental constraints from 
Ref. [6] - see the text for details. 

do not predict any mass loss (APR and SkM*) or they predict a mass loss 
which is too large. 

Another interesting possibility has been recently suggested by Lattimer and 
Schultz [6] . Measurement of the moment of inertia of Pulsar A (with known 
mass of about 1.34 M Q ) in the same system, to ~10% precision, would allow a 
rather accurate estimate of the radius of the star and of the pressure in matter 
of density in the vicinity of 1 -j- 2uq. Such information could provide a rather 
strong constraint on the EoS of neutron stars. In Fig. 7, using exactly the same 
scale as as Fig. 1 of Ref. [6], we illustrate this possibility. The dashed horizontal 
lines depict the error range of a hypothetical result of a measurement, selected 
as close as possible to the expected value. If this were the actual experimental 
result, many EoS would not be useful in determining the radius of Pulsar A, 
as they would cross the area between the two dashed lines either below or 
above M = 1.34 M Q . The QMC EoS would satisfy this constraint. Moreover, 
their relatively weak mass dependence would be an advantage in the study of 
of another analogous system with known mass up to ~1.6 M Q (N-QMC7r4) 
and ~2 M . We note that incorporation of hyperons shortens this range at 
higher masses, as illustrated in Fig. 7. 

Finally, we study the implications of the baryon composition of stellar mat- 
ter, as calculated using the sample EoS, for a possible cooling mechanism of 
neutron stars either just after their birth in supernovae or after heating in an 
accretion episode. The cooling mechanisms can, in principle, provide further 
important constraints on neutron-star models. However, as discussed in more 
detail previously [44], the cooling processes for both young and old neutron 
stars are not currently known with any certainty, although several scenarios 
have been proposed [49,50]. Many of them involve emission of neutrinos from 
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Fig. 8. The proton fraction y v for beta stable matter, as calculated using the QMC 
EoS is plotted as a function of the baryon number density ng. The data for the 
APR EoS are added for comparison. 

the stellar core. The most frequently discussed, the direct URCA process, pro- 
duces neutrinos in proton-neutron weak decays with an additional production 
of either electrons of muons. The proton concentration in the matter is a cru- 
cial ingredient of the process. Conservation of energy and momentum requires 
Vr/ 3 < y]J 3 +Z/y 3 ; or vl/ 3 < y^! 3 + yj/ 3 i to be satisfied. This, in turn, can hap- 
pen only in matter where the proton fraction grows steadily with increasing 
density and y p is greater than about 0.11 for matter consisting of only protons, 
neutrons and electrons. Medium effects and additional interactions among the 
particles change this number only slightly [46] but the presence of muons in- 
creases the critical ratio, y p , to about 0.15. We display the proton fraction for 
sample EoS in Fig. 8. It is clear that the proton fraction in F-QMC700 and 
F-QMC7r4 pass this critical limit at densities close to 0.6 fm -3 , coincident with 
the appearence of negative S~ hyperons. This density is well below the central 
density for maximum mass models and thus the direct URCA process is al- 
lowed for models of astrophysical interest, in agreement with expectation. On 
the contrary, for the N-QMC700 and N-QMC7r4 models the critical density for 
the start of the direct URCA process is much higher, between 0.8 -j- 1.0 fm -3 , 
and is close to or exceeds the inferred maximum central density, suggesting 
that this cooling mechanism would not be allowed in these cases. This re- 
inforces our message concerning the importance of incorporating of the full 
baryon octet in QMC. 



6 Conclusions 



One of the most important features of the QMC EoS is that although the QMC 
and Skyrme energy functionals have similar structure at densities close to 
nuclear saturation density, they differ significantly at higher densities, because 
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of their different density and isospin dependence. The QMC EoS provide cold 
neutron star models with maximum mass 1.9-2. f M , with central density less 
than 6 times nuclear saturation density and thus offer a consistent description 
of the stellar mass up to this density limit, without extrapolation beyond the 
region of validity - as needed, for example, when the non-relativistic, nucleon- 
only Skyrme models are used at these densities. The maximum mass of neutron 
stars in the QMC model is close to the most recent observation, the central 
density is relatively low in comparison to other models and all the predicted 
stars are well within the causal limit. 

The QMC model includes self-consistently the presence of hyperons at high 
densities. Unlike most other models, no hyperon contribution at densities lower 
than 3 no is predicted in /3-equilibrium matter. At higher densities, and 
A hyperons are present. The absence of the S ± '° hyperons is due, on the one 
hand, to the unique features of the model, namely the implementation of the 
color hyperfine interaction and the scalar polarisability of the baryons and, 
on the other hand, to the use of the Hartree-Fock approximation. The pro- 
nounced effects of the presence of hyperons on the properties of dense matter 
and neutron star models have been identified by a comparison of nucleon-only 
EoS with full-baryon-octet EoS. In nuclear matter, the presence of hyperons, 
as expected, softens the EoS as illustrated in Fig. 4. However, it needs to be 
stressed that the composition of stellar matter at high density is not neces- 
sarily the only determining factor of the behaviour of pressure with increasing 
density. In fact, it is critically dependent on the potentials acting amongst 
particles present in the matter. 

In neutron star models, the hyperon component systematically reduces the 
prediction of the maximum mass of a neutron star and increases its radius 
at maximum mass, as demonstrated in Fig. 5. This is a very important find- 
ing, because many non-relativistic EoS, including those based on the effective 
Skyrme interaction and the Al8+5v+UIX* potential, do not include hyperons 
at densities corresponding to the central density of maximum mass star models 
and thus their results should be taken with caution. Moreover, the relation- 
ship between the moment of inertia and the stellar gravitational mass, M g , is 
clearly influenced (see Fig. 7), as is the cooling mechanism expected to take 
place in neutron stars - the direct URCA process - is predicted to be allowed 
in stars with hyperons but not in nucleon-only neutron stars. It is also impor- 
tant to realise that the models for the nucleon-hyperon and hyperon-hyperon 
interactions play an important role in this comparison. In future, it would 
be interesting to apply the present QMC model to non-equilibrium matter 
at finite temperature and produce EoS that could be tested in core-collapse 
supernova simulation models. The high density threshold for the appearance 
of hyperons probably rules out any hyperon-related effects taking place in the 
collapsing star at the bounce density of 1.5 — 2.0 no, in variance with some 
other EoS with a much lower threshold density (e.g. [51]). 
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8 Appendix 

8. 1 Free mass of the baryons in the bag model 

The free mass of the baryon of flavor / with quark content (N u d, N s ) is written 

as 

Mj = W + NMm.) _ Zo + ) + 4^ 

R R 3 

where the quark mode Q(m) is determined by the boundary condition 
sinx 1 / sinx 



+ - cosx =0, x = JVL 2 - (mR) 2 (45) 



x fl + mR \ x 



and B is the bag constant. In the presence of the scalar field one can have 
Q 2 — (mR) 2 < but the equation remains valid by analytical continuation. 
The bag radius is fixed for each flavour by the stability condition 



dM f 
~~dR 



0. (46) 



We set the u, d mass equal to zero and the zero point parameter Zq is assumed 
to be the same for all particles. For completeness we briefly recall how the 
hyperfine color interaction AEm is evaluated according to [52p"l. For a color 
singlet baryon one has 

AE M = Y i E ij a i .<? j (47) 

i<j 



4 We thank H. Grigorian who pointed out a missprint in Equation 2.23 of Ref.[52] 
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Table 7 

Spin-flavor factors for Eq.(54) 

with Gi the Pauli spin matrices of the i th quark and 



E 



i.i 



IJt,i{ R)iij{R) 
R 3 



(48) 



where a c is the color coupling constant. The magnetic moment has the ex- 
pression 

= R 4flK) + 2i?m t - 3 
W ; 6 2n(mi) (fiK) - 1) + i?m/ 1 j 



and we point out that its quark mass dependence will produce a non trivial 
flavor dependence of the coupling to the nuclear scalar field. The expression 
for the overlap integral 1^ is: 



^ ^ ^UiUj Ax^Xj sin sin Xj -\- x^x j K%j (50) 
2(xi sin 2 Xi - lyi){xj sin 2 ^ -\yj) 



Kij = 2x i Si{2x i ) + 2x j Si(2x j ) (51) 

— (xi + Xj)Si{2xi + 2xj) — (xj — Xj)Si{2xi — 2xj) (52) 

Hi = Xi — sinxj cosxj, X; = -y/f2(mj) 2 — (Rrrii) 2 (53) 
Using the spin-flavor wave functions of the baryons one can write: 

AE M = a 00 E 00 + a 0s E 0s + a ss E ss (54) 



where the index refers to the u, d quarks ans s to the strange quark. The 
spin-flavor factors a„ are given in Table 7. 

To fix the free parameters (B, Z ,m s ,a c ) we take the free nucleon radius 
R^ ee as the free parameter and require that the nucleon and A(1232) have 
their physical massed"^! Together with the stability equations this determines 

5 The expression for the energy of A(1232) is the same as for the nucleon except 
that ago = 3. 
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R f N ee ( fm) 


B{ fm" 4 ) 


Zo 


a c 


m s ( MeV) 


0.8 


0.554 


2.642 


0.448 


341 


1.0 


0.284 


1.770 


0.560 


326 



Table 8 

Values of the parameters (B, Zo,m s ,a c ) for typical values of R^ ee ■ 

(B, Z , a c ). We then choose m s so as to have a best fit to the (A, S) masses. 
The resulting parameters are shown in Table 8 for typical values of R^ ee . 



8.2 Effective mass 



We now consider the effect of the nuclear scalar field which is assumed to 
be constant over the volume of the baryon. We recall that the variation of 
the field just produces the spin-orbit interaction, which we can neglect when 
considering uniform matter. Let a be the uniform value of the scalar field and 
gl its coupling to the u, d quarks. We assume that the strange quark does 
not interact with the scalar field. Therefore the coupling to the scalar field 
amounts to the substitution 



m ud 



m ud 



9>, 



m. 



(55) 



in the free mass equation, Eq. (44), followed by application of the stability 
condition to determine the actual radius. The mass then becomes a function 
of cr and of the parameters, (B, Z , m s , a c ), which themselves depend only on 
Rpj ee '■ So we have 



My (in medium) = Mf(m u d 



(56) 



Since we have two independent energy scales, a and 1/R^ ee , dimensional anal- 



? free- ) 
N 



ysis is not very helpful and we simply make a fit of Mj(a, Rj^ ee )—Mf(Q, R 
in powers of m U( i— g q a a ■, with the coefficients also fitted as polynomials in R^ ee . 
We then define the free o — N coupling as 



dM N (a,R{r 



9a 



N 



da 



-9l 



dM, 



N 



'a, R{ ree " 



dm 



(57) 



m=0 



which allows us to eliminate g% in favor of g 
the form 



We then get an expansion of 



M f (a,R f D 



M f (0,R f D 



P 



{1) (R f N ree )9*° 



P 



m {R f DM 2 



....(58) 
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where, by construction 



P {1) (R 



N 



-1. 



(59) 



If the mass were approximated by 
N u4 tt(m ud ) + N s tt(m s ) 



M f 



TDfree 



(60) 



we would have 

= -2/3, P| 1} 



-1/3, 



(61) 



but this approximation is severely broken in the realistic case. The order of the 
polynomials in R has been set to 2 and the coefficients have been fitted over 
the range [0.7 fm — > 1.3 fm]. We have checked that an expansion truncated at 
order (g a a) 2 was sufficient, even for scalar fields as large as g a a = 600 MeV. 
This corresponds to densities of order 1 fm -3 , which is sufficient for our pur- 
poses. The expansion becomes, with everything expressed in fm: 



M N (a) = M N - g a a 



0.0022 + 0.1055i$ ee - 0.0178 (R^ ee 



M A (a) = M A 



0.9957 - 0.22737i?^ ee + 0.01 (i$ ee ) 
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+ 



0.0022 + 0.1235i?{T e - 0.0415 (R s Z ee 



M A (a) = M A 



0.6672 + 0.0462i2jf" ee - 0.0021 (i^ ree ) 



2 



9aV 
2 



+ 



0.0016 + 0.0686i^ ree - 0.0084 (R% 



M s (a) = M s 



M*(a) = M* 



0.6706 - 0.0638i?^ ree - 0.008 (i?^ ee ) 



0.0007 + 0.0786i^ ree - 0.0181 {R^ ee 



2 



9a<y 

2 



0.3395 + 0.02822i^ ree - 0.0128 (i?^ ee ) 



9aO 

(9<r<r) 

2' 



0.0014 + 0.0416i$ ee - 0.0061 (R% ee 



1 iY 



9o<J 

2 



(62) 



(63) 



(64) 



(65) 



(66) 



Note that in these calculations we use the physical masses for the free par- 
ticles. The influence of the bag model dependence then appears only in the 
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N 


A 


s 




d( fm) 


0.15 


0.15 


0.15 


0.15 


w f 


1 


0.703 


0.614 


0.353 




1 


0.68 


0.673 


0.371 



Table 9 

The weights for R^ ee = 0.8 fm. 

interaction piece, Mf(a) — Mj(0). It is convenient to write the dynamical mass 
in the form 

M f (a) = Mf — w a f g a a + d -w) (g a af (67) 



where the scalar polarisability, d, and the dimensionless weights, wj,wj, are 
deduced from Eqs. (62-66). For illustrative purposes we give their values in 
Table 9 for our preferred value of the free radius, R^ ee = 0.8 fm. 
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